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We propose the Wilson discretization of the supersymmetric Yang-Mills Quantum Mechanics as a lattice version 
of the matrix model of M-theory. An SU(2) model is studied numerically in the quenched approximation for D=4. 
A clear signal for the existence of two different phases is found and the continuum pseudocritical temperature is 
determined. We have also extracted the continuum limit of the total size of the system in both phases and for 
different temperatures. 



1. Introduction. Since the discovery of the web of 
duahties linking different superstring theories, it 
has been suggested that all these theories are just 
different perturbative expansions of a single II- 
dimensional theory christened M-theory, The true 
nature of this theory remains, however, mysteri- 
ous. In Banks, Fischler, Susskind and Shenker 
proposed a concrete formulation for the degrees 
of freedom of M-theory in a flat Il-dimensional 
background. 

In this contribution we would like to report on 
the work done in on a lattice formulation of the 
BFSS matrix model, thus allowing for a nonper- 
turbative treatment. We use the version with eu- 
clidean signature suitable for studying thermody- 
namical properties of the theory. In fact it turns 
out that there is a very rich spectrum of physics 
involved. A study of these properties forms the 
main motivation for our investigation. Below we 
briefly review the physics involved as it also pro- 
vides an incentive for further study. 

The BFSS lagrangian for flnite iV is a dimen- 
sional reduction of lOD SU(N) SYM down to I 
dimension, schematically 

S^j dt (^^{DtX'f + [X^X^]"^ + fermionsyil) 

Equivalently it describes a system of N DO 
branes in the decoupling limit. All the dynam- 
ics at finite temperature T can be expressed in 
terms of a dimensionless coupling constant g^^j = 
gyj^N /T^ . At high temperatures the system is 
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perturbative. As we lower the temperature we 
enter the strong coupling regime. 

At this stage we may apply the framework of 
the AdS/CFT correspondence (see [||,|| in this 
context). An equivalent description is provided 
by a 10 dimensional Black Hole (BH) supergrav- 
ity solution, whose Hawking temperature is iden- 
tified with the temperature of the SYM quantum 
mechanics (^. In particular the Bekenstein en- 
tropy of the black hole calculated from the area 
of the horizon should be identified with the en- 
tropy of the SYM quantum mechanics (for recent 
new results see |^ ) which thus gives a microscopic 
realization of the 'macroscopic' BH entropy. It 
would be very interesting to directly observe this 
prediction by studying numerically (|l]). 

Another interesting point is the nature of the 
transition from the perturbative to the strong 
coupling phase. In \^ it has been identified with 
the Horowitz-Polchinsky correspondence point 
between a black hole phase and 'excited string' 
phase. It's detailed nature is still unknown. 

As one lowers the temperature, one reaches a 
regime where the string coupling becomes strong 
and one has to lift the lOD BH solution to a IID 
'black wave'. Further down on the temperature 
scale the IID 'black wave' localizes to a 1 ID BH. 

All the above physics should be contained in 
the strong coupled regime of (0). For most of 
these questions supersymmetry is crucial. In 
our exploratory study we perform quenched sim- 
ulations and focus on showing the existence of 
two phases in the high and low temperature 
regimes. When performing simulations in the full 
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unquenched model the observed phases should 
correspond to the black hole phase and the per- 
turbative phase. In the following we present the 
lattice formulation, give details on the algorithms 
and present results of first exploratory quenched 
simulations for = 2. 

2. Lattice formulation. A direct Monte-Carlo 
simulation of the action is rather nontriv- 
ial. The potential for the bosonic fields [X^,X^]'^ 
has flat valleys which may be difficult to handle 
numerically. Furthermore it is essential to im- 
pose the Gauss law constraint or equivalently to 
project onto gauge invariant states. This may be 
important especially as the supersymmetry alge- 
bra in the continuum closes only on the space of 
gauge invariant states. Nevertheless the nonper- 
turbative studies of the IIB matrix models (SYM 
reduced to zero dimensions) in this formulation 
are well advanced ||J^. 

In order to overcome the above difficulties, and 
to use standard, well developed lattice techniques 
we decided, however, to follow a different route. 
The action comes from the standard SYM ac- 
tion in D = 10, upon assuming that all fields 
do not depend on spatial coordinates. We may 
now use the Wilsonian discretization in D = 10 
and impose the constancy of link variables in 
the spatial directions. In this manner we ob- 
tain at once a D dimensional hypercubic lattice 
A^i X . . . X Nd reduced in all space directions to 
Ni — l,i — 1...D — 1. Gauge and fermionic vari- 
ables are assigned to links and sites of the new 
elongated lattice in the standard manner. The 
gauge part of the action is 

= E E ^^e(TrC/^.(m)), (2) 

with /? = 2N/a^g^, and Uf,^{m) = Ul{m)Ul{m + 
u)U^(m + fi)U^{m), U^{m) = exp {iagA^{am)), 
where a denotes the lattice constant and g is 
the gauge coupling in one dimension. The in- 
teger time coordinate along the lattice is ni. Pe- 
riodic boundary conditions Uf_i {m + i^) = C/^ (to) , 
= 1 . . . Z? — 1 , guarantee correct classical contin- 
uum limit (|l]) with Xi = gAi. In this formulation 
the projection on gauge invariant states is natu- 
rally implemented. Addjoint fcrmions can be in- 
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fit: 


/3c = 





X^/NDF a 7 
■55/2 .17 ±.05 3.02 ±.33 



Table 1 

Estimated location of the transition region f3c S 
{Plow, Pup) for different lattice sizes Nt and results 
of the power fit. 



eluded as in the lattice studies of SUSYM theory 
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2. Results. We choose the distribution of the 
Polyakov line P = -^TrJI^'^^ Uuim), as an or- 
der parameter. Similarly to lattice QCD, sym- 
metric concentration of the eigenvalues around 
indicates a low temperature phase (which would 
have here the interpretation of a black hole phase) 
where (P) ~ 0, while clustering around ±1 (for 
SU(2)) is characteristic of the high temperature 
(elementary excitations) phase. Our first result 
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Figure 1. Distribution of the Polyakov line P, 
-1 < P < 1, for different /3 and Nt. 



is evident from Fig. 1. For each Nt a definite 
change of the shape with /3 is observed. The sys- 
tem (^) definitely shows the onset of the phase 
change, even in the quenched approximation and 
for N = 2. Second, the dependence of the pseud- 
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ocritical temperature Pc on the time extent Nt is 
consistent with the continuum hmit expectations 
Tc ~ (g^TV)^/^ !§. Indeed, the temperature of 
a system is given by T = l/{aNt) which impHes 
Pc ~ Nf. This is confirmed in Table I where jS in- 
tervals where the transition occurs are presented 
for several lattice sizes Nt . Results of the power 
law fit are also quoted. A good quality of the fit 
and the agreement with the canonical exponent, 
7 = 3, is encouraging. Simultaneously, we obtain 
the proportionality coefficient a which gives 

- i^)'^'i9'Ny/' = (.28±.03)(ff2iV)i/3.(3) 

The coefficient in this relation has been deter- 
mined for the first time. Only proportionality of 
the two scales was considered [||J^,^- The pseu- 
docritical temperature and a may depend on N, 
therefore similar analysis for higher gauge groups 
is necessary. 

Next we study the total size of the system — 
5^Ea(^")^ ®- We define for SU{2) 

(i^2)^(4-((Tr(C/.)^)/a^ (4) 

where Us is any space link. 

We used mostly the standard local Metropo- 
lis update with enough thermalization and decor- 
relation sweeps to take care of the critical slow- 
ing down. For example when running at d — 1.0 
we used 5000 thermalization and 50 decorrelation 
sweeps, while for d = ag^^^ — 0.1 about 10® ther- 
malization and 5000 decorrelation sweeps were 
required to achieve independence of the starting 
configuration. This also agrees with the dynam- 
ical exponent z = 2. In addition we have devel- 
oped the new heat bath algorithm designed for an 
update of the space-space plaquettes in (|^) which 
contain twice the same SU(2) link. In the stan- 
dard S3 parametrization of SU(2) distribution of 
the three vector u,\u\ — u < 1 is the product 
of the three dimensional gaussian and the cosh 
factor which depends on y^l — u^). At fixed val- 
ues of the neighbouring links the center of the 
gaussian may lie outside the kinematical region 
u < 1, which even for intermediate /3 makes the 
simple gaussian generation inefficient. However 
by a suitable choice of the coordinate system^ 

■^The one with the third axis pointing towards the maxi- 



the problem becomes manageable. The quadratic 
part of the log cosh factor is included in the above 
gaussian and deviations from the exact function 
are included by the accept-reject step. Results 
obtained with the new heat bath and the stan- 
dard Metropolis agree within statistical errors for 
(3 < 64. For higher /? the accept-reject step be- 
comes inefficient. For higher /3 we have also mon- 
itored the correlation length in the torelon chan- 
nel at zero temperature. It shows the canonical 
scaling with a as expected. Fig. 2 shows the de- 
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Figure 2. Dependence of the total size of the 
system on a, for f = Tg''^'^ = .1 - .3, .6, .9, 1.2 
and 1.5 (upwards) in units of g^^^- Quartic fits 
are represented by the solid lines. 



pendence of < i?^ > (in units of 5^^^) on a, for 
several values of the temperature T. MC results 
depend smoothly on a, at fixed T, confirming the 
existence of the continuum limit (^. The a de- 
pendence is different in low and high temperature 
regions. For T > 1.5 simulations for smaller a are 
required in order to see the quadratic approach to 
the continuum. Practically the same results were 
obtained when we extracted {J^) from another 
lattice observable \Tr{Us)\- Fig. 3 shows the size 
of a system extrapolated to a = as a function 
of the temperature. Both quadratic and quar- 
tic fits of a dependence were used. The stability 

mum of a gaussian on a boundary sphere u = 1. 
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Figure 3. Size of the system extrapolated to 
the continuum, as a function of the temperature. 



of quadratic fits with respect to removing one or 
two data points with smallest a was also checked. 
Results of the extrapolation were stable with re- 
spect to all these variations. Displayed errors in- 
clude above systematic effects. The location of 
the transition region agrees approximately with 
the estimate (^) of the pseudocritical tempera- 
ture Tc = 0.35 ± 0.04 Again, it is evident that 
the system is indeed different in the two regimes. 
Moreover, our results agree qualitatively with the 
analytical prediction obtained from a gap equa- 
tion in the infinite N limit ||]. The latter gives 
a temperature independent constant at low tem- 
peratures and the classical T^/"^ growth for high 
temperatures. We have also found a reasonable 
agreement with a simple mean field model for 
SU{2) with the gauge projection which will be 
discussed elsewhere. 

3. Outlook. We have constructed, for the first 
time, a matrix model of M-theory on a lattice. 
MC simulations show that even the strongly sim- 
plified version of the model has different high and 
low temperature phases as expected in the orig- 
inal theory. Pseudocritical temperature and to- 
tal size of the system show canonical approach 
to the continuum limit. This indicates that, at 
least in this respect, the system may be simpler 
than the ones encountered in lattice studies of 
field theories in extended space-time. In particu- 



lar, restoration of the supersymmetry, broken by 
lattice discretization, is feasible. 

Of course results presented here are only a hint 
that we may be on a right track. One urgent task 
is to repeat present programme for higher N and 
D=10. This, although relatively simple, may re- 
quire some refinement of existing algorithms or 
inventing new ones dedicated to the linear sys- 
tems. Equally important is to include dynamical, 
supersymmeric fermions. For D=4 this is again 
relatively standard and would allow us to study 
full intricacies of lattice breaking and subsequent 
restoration of the supersymmetry in the contin- 
uum limit. For D=10 the problem of dynamical 
fermions is open and presents an exciting chal- 
lenge. Again the linear nature of the system may 
offer some simplifications. 

To conclude, we believe that many interesting 
applications became available opening new pos- 
sibilities for the nonperturbative studies of the 
prototypes of M-theory. 
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